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The structure of the (v5 x v / 5)i?26.6° reconstruction of LaAlC>3 (001) has been determined using trans- 
mission electron diffraction combined with direct methods. The structure is relatively simple, consisting of a 
lanthanum oxide termination with one lanthanum cation vacancy per surface unit cell. The electronic structure 
is unusual since a fractional number of holes or atomic occupancies per surface unit cell are required to achieve 
charge neutrality. Density functional calculations indicate that the charge compensation mechanism occurs by 
means of highly delocalized holes. The surface contains no oxygen vacancies and with a better than 99% confi- 
dence level, the holes are not filled with hydrogen. The reconstruction can be understood in terms of expulsion 
of the more electropositive cation from the surface followed by an increased covalency between the remaining 
surface lanthanum atoms and adjacent oxygen atoms. 

PACS numbers: 68.35.Bs, 61.14.Lj, 68.37.Lp, 73.25.+i 



Lanthanum aluminate (LaA103, LAO) is an important ox- 
ide material as a substrate for thin film growth |lL [A [j, |4J] a 
potential gate dielectric (or buffer layer) Bi |rjj U H Eft UoQllll , 
as well as a catalyst This transition metal oxide 

is representative of a larger class of materials with ABO3 sto- 
ichiometry known as perovskites. There has recently been in- 
terest in the character of the interfaces between LaAlC<3 (001) 
and other materials such as Si and SrTi0 3 E [H. 

Since 

both La and Al are almost always in a trivalent state, conven- 
tional methods of counting carriers for charge compensation 
lead to the unusual conclusion that a fractional number of ex- 
cess carriers per unit cell is required. The effect of such cor- 
related electron properties in oxide heterostructures has been 
studied extensively theoretically Il7lll8lll9ll20ll21ll for these 
interfaces, since charge redistribution is well-described by the 
electrochemical potential gradient between the materials. 

LaA103 (001) consists of alternating layers of LaO and 
AIO2 stacked along the <001> cubic direction. Conse- 
quently, the formal charges of La 3+ , Al 3+ and (D 2 ~ pro- 
duce two terminations differing in nominal charges of (La- 
0) + and (Al-C<2)~, and an excess half electron (or hole) 
exists per unit interface cell. The LAO (001) surface is 
therefore polar and classified as a Type III surface within 
Tasker's convention [22]. The composition, structure and 
morphology of LAO (001) has been researched in the past 
few years; however, structural information of how the surface 
terminates is still ambiguous. The majority of the analyses 
flmHmS^SllilB have considered a (1 x 1) surface 
prepared by relatively low-temperature annealing (the melt- 
ing point of LaA103 is 2110°C) in various atmospheres, and 
report either (or both) LaO and AIO2 terminations; it is not 
clear whether the surfaces have reached thermodynamic equi- 
librium. The only reported reconstructions on the (001) sur- 
face of LAO is a (5x5) reconstruction obtained by annealing 
at 1500°C for 20 hours in air Q . 

To date, all models have failed to offer any feasible sta- 
bilization pathways that are consistent with the experimental 
results. Furthermore, very little of the work has addressed the 



issue of charge compensation and how to reconcile the dipole 
at the surface. In this Letter, we report experimental results 
for the {V5xV5)R26.6° surface reconstruction, hereafter re- 
ferred to as Rt5, on the LaAlO3(001) surface. We also pro- 
vide a first-principles investigation of the structure and show 
how the polar surface is passivated through a redistribution of 
the near surface electron density. 

Single crystal LaA103 (001) wafers from MTI Corporation 
(99.95% pure) were prepared for transmission electron mi- 
croscopy studies using conventional methods of dimpling and 
ion-beam thinning (a Gatan Precision Ion Polishing System 
with 4.8 kV Ar + ions) until electron transparent. They were 
then annealed at temperatures between 1 100-1500°C in a Car- 
bolite STF 15/51/180 tube furnace for three hours. While ini- 
tial experiments were performed in air, we also annealed in a 
mix of 20% O2 : 80% N2, which mass-spectrometer measure- 
ments indicated had a maximum impurity level of 10 ppb of 
H 2 0. 

Diffraction experiments were performed with a Hitachi 
8 100 electron microscope operating at 200 kV. Negatives with 
exposure times ranging from 0.5 to 90 s were recorded then 
digitized to 8 bits with a 25 jim pixel size on an Optron- 
ics P-1000 microdensitometer which was calibrated to be lin- 
ear over the selected exposure range. Intensities measure- 
ments were determined using a cross-correlation technique 
03111 . and the data sets (7911 surface reflections) were sym- 
metry reduced using p4 plane group symmetry to 94 indepen- 
dent beams with the error for each reflection determined using 
robust statistical methods. These were analyzed using elec- 
tron direct methods (EDM 2.0.1) software 113211 and in-plane 
atomic positions refined against the experimental data gave a 
r 4.83. 

For the theoretical modeling, the LAO surface structure was 
geometry optimized using a three-dimensional periodic DFT 
surface slab model of 9 layers (118 atoms) separated by 8 A of 
vacuum. To analyze the charge density, calculations were 
performed using the all-electron (linearized) augmented-plane 
wave + local orbitals (L/APW+lo) method as implemented in 





FIG. 1: (color online). Dark field image showing extended <100> 
faceting with step bunches and reconstructed terraces. A small probe 
off-zone diffraction pattern (inset) is shown with the rhombohedral 
(lxl) bulk unit cell (orange, right) and the surface unit cell for both 
domains of the reconstruction (green, left). 




FIG. 2: (color online). Top view of the Rt5 surface reconstruction 
(left) with the surface unit cell outlined in black (a=8.537 A). Con- 
tour map of the valence density in the surface plane (right) at 0. 1 e/A 3 
showing the increased electron density between the surface oxygen 
(01 and 02) and lanthanum atoms (Lai) corresponding to increased 
covalent bonding. Atoms 01, 02 and Lai are in the surface layer, 
while All, 03, 04 and 05 are in the second layer. 



WIEN2k Q with the PBE-GGA Q exchange-correlation 
functional, a plane-wave cutoff of RK max =6.75 and muffin-tin 
radii of 1.75, 1,75 and 2.33 Bohr for O, Al and La respec- 
tively. To test for water splitting and oxygen vacancies, cal- 
culations were also performed using the projector augmented 
wave (PAW) approach ll35Tl as implemented in the Vienna 
Ab Initio Simulation Package (VASP) code |36[ 37] using a 
3x3x1 k-point grid, plane wave energy cutoff of 360 e V, elec- 
tronic iteration convergence of 0.001 eV and geometry relax- 
ation convergence of 0.01 eV. The Fermi surface was smeared 
using a Gaussian width of 0.20 eV. For these calculations we 
used the conventional LDA energy functional as well as the 
PBE Q and PW92 Q functionals. In all cases the surface 
unit cell used was matched to the DFT minimized lattice pa- 
rameter for the appropriate functionals. 

A typical transmission electron diffraction pattern of the re- 
constructed surface is shown in Fig. [T] The surface is well 
ordered with facets evident in the dark field image and mini- 
mal diffuse background scattering from disorder present in the 
diffraction pattern (inset). Direct methods analysis and refine- 
ment were straightforward, and gave the structure shown in 
Fig. |2 information about the atomic positions can be found in 
Table H] while the experimental and calculated intensities are 
given in the supplemental material. 

Bulk LAO has a rhombohedral to cubic phase transition at 
435±25°C. Consequently, the Rt5 reconstruction forms on a 
cubic LAO substrate, but this substrate becomes rhombohe- 
dral upon cooling. As the structures of rhombohedral LAO 
and cubic LAO only differ by tenths of a degree, it is rea- 
sonable to expect the cubic surface reconstruction to be ac- 
commodated on bulk rhombohedral LAO at room tempera- 
ture. Nonetheless, the possibility of a rhombohedral surface 
structure (p2 plane symmetry) and a cubic surface structure 
(p4, p4mm, p4mg plane symmetry) were both explored and a 
p4 symmetry provided the best fit to the data. The possibility 



of an aluminum overlayer was also considered, however the 
fit to the data was noticeably inferior. 

The Rt5 reconstruction is an overlayer of lanthanum oxide 
on the AIO2 bulk termination of LAO with one lanthanum va- 
cancy per surface unit cell. The surface stoichiometry may be 
written as (VLa^s) - ?, where V is the lanthanum cation va- 
cancy and the surface unit has a nominal charge of — A. Each 
surface lanthanum (Lai ) is coordinated to four oxygens within 
the surface layer and four oxygens in the layer below. Ad- 
ditionally, the lanthanum atom is displaced into the bulk by 
«0.20A, while the oxygen atoms are displaced away from the 
bulk. The oxygen atoms (Ol and 02) in the surface layer pre- 
serve the octahedral coordination of the aluminum atoms in 
the layer below; however due to the La vacancy, Ol is now 
only four-coordinate while 02 is five-coordinate. 

The formation of the Rt5 reconstruction can be understood 
as follows. To reduce positive charge on the surface of an ideal 
La-O bulk termination, either La 3+ or Al 3+ cation vacancies 
are formed. Because La 3+ is less electronegative than Al 3+ , 
lanthaum cation vacancies are favored. It follows that the 
surface bond covalency should increase to satisfy the under- 
coordination of the surface oxygens. Therefore since the La-O 
bonds are longer compared to the Al-O bonds, and the non- 
bonding oxygen repulsive interactions are smaller than those 
of the Al-O octahedra, the La-O bonds become shorter, with 
an average La-O bond length of 2.61 A (compared to 2.68 A in 
the bulk). However, despite these significant changes, a per- 
sistent, albeit reduced, polarization remains at the surface. 

The lanthanum cation expulsion nearly reconciles the 
charge neutrality problem at the surface. For the Rt5 re- 
construction with an area equal to five interface unit cells 
of the bulk, it is impossible to form a fully charge compen- 
sated surface without invoking a fractional density (per re- 
constructed surface cell) of carriers or a partial occupancy of 
cation sites. While there was too little diffuse intensity in the 
diffraction data to support a fractional occupancy of sites, a 
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TABLE I: Atom positions for the Rt5 reconstruction are given as 
fractional coordinates of the surface cell (a = 8.526 A). The La va- 
cancy is denoted by V at the surface (layer 1) and the bulk corre- 
sponds to layer 5. A z = |zdft ~ ZBulkl m A and positive (negative) 
deviations indicate a displacement away from (into) the bulk. The 
excess charge values for each atom are denoted by <5Q (10~ 2 e), and 
the integrated hole density as p(h) (10~ 2 e/A 2 ). Bulk LAO provides 
the charge reference: O = — 1.540e, La = 2.065e and Al = 2.554e. 



Experiment 


WIEN2K 






Layer Atom x y 


x y z 


A 2 


5Q p(h) 



1 


V 


0.00 


0.00 


0.000 


0.000 


0.320 


— 


— 


— 




01 


0.34 


0.90 


0.344 


0.898 


0.320 


+0.036 


10.2 


2.1 




02 


0.50 


0.50 


0.500 


0.500 


0.321 


+0.055 


5.8 


2.5 




Lai 


0.81 


0.40 


0.808 


0.396 


0.310 


-0.196 


-2.6 


0.3 


2 


All 


0.30 


0.90 


0.298 


0.901 


0.240 


+0.035 


-0.7 


0.0 




A12 


0.50 


0.50 


0.500 


0.500 


0.239 


+0.015 


-0.7 


0.0 




03 


0.50 


0.00 


0.500 


0.000 


0.232 


-0.169 


4.1 


2.0 




04 


0.21 


0.11 


0.210 


0.106 


0.239 


+0.009 


1.4 


0.6 




05 


0.70 


0.59 


0.702 


0.594 


0.237 


-0.041 


2.0 


1.3 


3 


06 


0.29 


0.90 


0.290 


0.899 


0.159 


-0.007 


2.0 


1.8 




07 


0.50 


0.50 


0.500 


0.500 


0.158 


-0.018 


2.2 


2.5 




La2 


0.80 


0.40 


0.800 


0.399 


0.158 


-0.036 


-0.3 


0.2 




La3 


0.00 


0.00 


0.000 


0.000 


0.163 


+0.099 


1.9 


0.2 


4 


A13 


0.30 


0.90 


0.298 


0.900 


0.079 


-0.003 


0.0 


0.0 




A14 


0.50 


0.50 


0.500 


0.500 


0.079 


-0.008 


0.3 


0.0 




08 


0.50 


0.00 


0.500 


0.000 


0.081 


+0.038 


0.7 


0.7 




09 


0.20 


0.10 


0.201 


0.100 


0.078 


-0.045 


0.4 


1.3 




010 


0.70 


0.60 


0.700 


0.600 


0.079 


-0.009 


0.5 


1.0 


5 


Oil 


0.29 


0.90 


0.304 


0.900 


0.000 


n/a 


1.3 


1.9 




012 


0.50 


0.50 


0.500 


0.500 


0.000 


n/a 


0.4 


1.0 




La4 


0.80 


0.40 


0.800 


0.400 


0.000 


n/a 


-0.2 


0.2 




La5 


0.00 


0.00 


0.000 


0.000 


0.000 


n/a 


0.3 


0.2 



low density of oxygen vacancies is undetectable experimen- 
tally. Similarly there is the possibility of disordered hydrogen 
atoms on the surface and several studies indicate that some 
bulk oxides (and surfaces) can contain low levels of hydrogen 

(S El SIH Si H). 

To test for the presence of hydrogen on the surface, we con- 
sidered the possible reaction 



4 (Rt5) + H 2 +^ 4 (Rt5H 



1/2 J 



?02 



with partial occupancy of a hydrogen bonded to Ol (lower in 
energy than the alternative 02). Four calculations were per- 
formed using VASP: one without hydrogen for a \/T0 x ylO 
supercell (9-layers, 236 atoms) rotated by 45° containing two 
surfaces each with four Rt5 unit cells, another of the same 
cell but with I hydrogen per Rt5 surface unit cell (atomic 
positions are given in the supplemental material) and two for 
the isolated molecules H2O and O2. The DFT calculations in- 
dicated that this reaction can occur with an energy change of 
-0.99 eV at T=0 K for the PBE functional, -1.04 eV with the 
PW92 functional and -1.46 eV for the LDA. By comparing 
the PBE and PW92 numbers we can estimate that the intrinsic 
surface uncertainty B46I1 is small, about 0.15 eV. A reason- 
able estimate of the error in the energies is |£lda — ^pbe|/2, 
or 0.24 eV. Thus, at a 99% confidence level (3(7) the maxi- 
mum energy for this exothermic reaction would be —1.71 eV 
at T=0 K. Using tabulated values for the free-energy of water 
and oxygen [47] the reaction becomes endothermic and would 
require more than 3.64 eV to take place at the Rt5 formation 
temperature of 1200°C in an air atmosphere with less than 




FIG. 3: (color online). Integrated hole density plot along the [100] 
direction at 0.01 e/A 2 shown next to the crystal structure. 



10 ppb water impurity. Hence, to a better than 99% confidence 
level when the reconstruction is formed at high temperature, it 
is energetically unfavorable for the surface to split water and 
incorporate H + . On cooling and exposure to environmental 
humidity it is quite possible that the surface splits water ll48ll . 

Similarly, to test for the presence of oxygen vacancies on 
the Rt5 surface, we considered two reduced structures with a 
4 oxygen vacancy per unit cell at sites Ol and 02 (the struc- 
tures are available in the supplemental material). It was found 
that at T=0 K the lower energy reduced structure (vacancy at 
Ol) is 2.58 eV higher in energy than the Rt5 structure as pre- 
sented. At the reconstruction formation temperature, the oxy- 
gen chemical potential is large enough to reduce the surface 
with an energy gain of -0.29 eV; however within the error as- 
sociated with these computational methods, we posit the prob- 
ability for this reduction is small. Furthermore, as the sample 
is cooled the oxygen vacancies (if any) would be filled and the 
Rt5 surface is preserved. Again, we find that an alternative 
charge compensation mechanism is required, and suggest that 
an electron hole fulfills this requirement. 

Hole densities were calculated by integrating over the 
Bader volumes 149m 1 . Somewhat unusual is that the hole densi- 
ties do not decay off into the bulk, but rather the fractional hole 
is very delocalized over all the oxygen atoms, as indicated by 
the projected hole density of states along the [100] direction 
normal to the surface shown in Fig. [3] This behavior is atyp- 
ical for most bulk oxides since hole trapping occurs through 
acceptor defect sites in open lattice perovskite structures 15011 . 
This effect leads to small hole mobilities in bulk materials, 
and is reminiscent of the Mott-Wannier exciton model in more 
ionic crystals (despite the increased covalency here). 

Although the charge on the atoms in a solid cannot be 
uniquely partitioned, there are several theoretical models 
which allow for its estimation; we used Bader's atom-in- 
molecule (AIM) method EH with the WIEN2k densities. Ta- 
ble |T] includes the excess charges (5Q) calculated for the var- 
ious layers. There is a noticeable decrease in the charge on 
the surface oxygen atoms (Ol and 02) as well as the sur- 
face lanthanum (Lai), with smaller variations decaying off 
more rapidly into the bulk for the other oxygen atoms. In 
bulk LaAlOi the Al-O bond has more covalent character than 
the La-O bond; the charge densities at the (3,-1) bond criti- 
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cal points, which are typically used as a qualitative estimate, 
are sa 0.45 and 0.24 eV/A 3 for the two. At the surface there 
is essentially no change in the covalent character of the Al-O 
bonds, but the density between the La-O at the surface in- 
creases dramatically to a value of « 0.50 eV/A 3 for the clos- 
est oxygen atom which is only 2.37 A away compared to 
2.70 A in the bulk (with the DFT optimized lattice parame- 
ters). This suggests that some of the charge on the oxygen 
ions in the bulk transfers into the covalent bonds near the sur- 
face, and rather independent of this, it occurs in the presence 
of a highly delocalized hole. 

For the Rt5 reconstruction on the LaA103 (001) surface, 
the surface polarity is quenched by the presence of the elec- 
tron hole. This point has strong implications on the growth of 
nano-layer composites and the dynamics describing surface 
reconstructions, that is, interface passivation is highly sensi- 
tive to the experimental preparation. Furthermore, it is likely 
that oxygen vacancies can strongly influence the formation of 
charge carriers at such interfaces; while the results suggest 
that they are filled upon annealing, it is important to recog- 
nize that their electron donor character may behave differently 
at interfaces. We find that the driving force for reconstruction 
is the polar discontinuity at the surface and charge neutrality 
requires either (or both) atomic or electronic reconfiguration: 
the expulsion of the lanthanum cation and delocalized electron 
hole meets this challenge. 
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